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ABSTRACT 

We  report  the  successful  fitting  of  a  Roche  model,  with  a  surface  temperature  gradient  following  the  von  Zeipel 
gravity  darkening  law,  to  observations  of  Altair  made  with  the  N avy  Prototype  Optical  Interferometer.  We  confirm  the 
claim  by  Ohishi  et  al.  that  Altair  displays  an  asymmetric  intensity  distribution  due  to  rotation,  the  first  such  detection 
in  an  isolated  star.  Instrumental  effects  due  to  the  high  visible  flux  of  this  first  magnitude  star  appear  to  be  the  limiting 
factor  in  the  accuracy  of  this  fit,  which  nevertheless  indicates  that  Altair  is  rotating  at  0.90  ±  0.02  of  its  breakup 
(angular)  velocity.  Our  results  are  consistent  with  the  apparent  oblateness  found  by  van  Belle  et  al.  and  show  that  the 
true  oblateness  is  significantly  larger  owing  to  an  inclination  of  the  rotational  axis  of  ^64°  to  the  line  of  sight.  Of 
particular  interest,  we  conclude  that  instead  of  being  substantially  evolved  as  indicated  by  its  classification,  A7 IV- V, 
Altair  is  only  barely  off  the  zero-age  main  sequence  and  represents  a  good  example  of  the  difficulties  rotation  can 
introduce  in  the  interpretation  of  this  part  of  the  HR  diagram. 

Subject  headings:  stars:  imaging  —  stars:  individual  (a  Aql)  —  stars:  rotation  —  techniques:  interferometric 


1.  INTRODUCTION 

Altair  (variously  a  Aql,  53  Aql,  HR  7557,  HD  1 87642,  of  spec¬ 
tral  type  A7 IV-V)  is  one  of  the  brightest  stars  in  the  northern  sky, 
sharing  membership  in  the  “Summer  Triangle”  with  two  other 
notable  A  stars.  Unlike  Vega  and  Deneb,  Altair  shows  a  rather 
diffuse  spectmm  that  was  early  recognized  to  be  due  to  a  large 
projected  rotational  velocity  variously  estimated  at  242  km  s_1 
(Uesugi  &  Fukuda  1982), 8  217  km  s_1  (Royer  et  al.  2002),  and 
200  km  s_1  (Abt  &  Morrel  1995).  These  estimates  of  its  pro¬ 
jected  velocity,  a  lower  limit  for  the  tme  rotational  velocity,  are 
already  a  significant  fraction  of  the  breakup  velocity,  estimated 
near  400  km  s 

Altair  has  become  a  significant  object  for  understanding  the 
atmospheres  of  main-sequence  stars  at  masses  near  but  above 
that  of  the  Sun.  Specifically,  Altair  and  a  Cep  are  the  two  hottest 
stars  showing  Ly a  and  C  n  emission,  taken  as  indicators  of  a 
chromosphere  (Simon  et  al.  1994;  Walter  et  al.  1995).  The  ab¬ 
sence  of  these  indicators  at  earlier  spectral  types  is  taken  to  mean 
that  significant  convection  disappears  at  this  point  on  the  upper 
main  sequence.  We  note  that  a  Cep  also  has  a  high  projected 
rotation  velocity  with  estimates  of 246  km  s_1  (Uesugi  &  Fukuda 
1982),  196  km  s_1  (Royer  et  al.  2002),  and  180  km  s~!  (Abt  & 
Morrel  1995)  listed. 

Altair’s  known  high  rotation  rate  has  prompted  attempts  to 
measure  the  geometrical  effects  of  its  rotation  over  the  years, 
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starting  with  the  Intensity  Interferometer  (Hanbury  Brown  et  al. 
1974).  However,  it  was  not  until  the  near-lR  observations  with 
the  Palomar  Testbed  Interferometer  (PTI;  Colavita  et  al.  1999) 
by  van  Belle  et  al.  (2001)  that  a  significant  flattening  was  de¬ 
tected.  Comparison  to  classical  (von  Zeipel  1924)  Roche  models 
showed  the  flattening  was  completely  consistent  with  the  ob¬ 
served  projected  rotation. 

Although  this  agreement  between  theory  and  observation  is 
nothing  short  of  epochal,  it  is  incomplete.  Except  near  breakup, 
apparent  oblateness,  interpreted  through  Roche  theory,  displays 
the  same  degeneracy  between  equatorial  velocity  and  tilt  (inclina¬ 
tion)  as  the  apparent  rotation  velocity:  one  determines  the  quantity 
peq  sin  i  well,  but  not  the  two  separately.  Nor  can  one  determine 
the  sense  of  rotation,  pro-  or  retrograde.  Besides  providing  a  test 
of  the  flattening  predicted  by  theory,  oblateness  measurements 
do  yield  the  position  angle  of  the  angular  momentum  vector,  the 
projection  of  that  vector  on  the  plane  of  the  sky. 

In  addition  to  flattening,  von  Zeipel  (1924)  predicted  that  for 
moderate  rotation,  stellar  disks  would  display  variable  surface 
temperatures,  hotter  on  the  rotational  axes  and  cooler  at  the  equa¬ 
tor.  Specifically,  if  one  defines  a  local  effective  gravity  accounting 
for  centrifugal  acceleration,  then  the  local  effective  temperature  is 
related  to  the  effective  gravity  as  ex  geS ,  which  is  referred  to  as 
“gravity  darkening.”  With  sufficient  rotation  and  at  intermediate 
inclinations,  gravity  darkening  predicts  that  stellar  disks  will  dis¬ 
play  asymmetric  intensity  distributions. 

As  we  describe  below,  this  prediction  is  of  great  interest  in 
the  field  of  optical  interferometry.  Asymmetric  intensity  distri¬ 
butions  produce  significant  imaginary  components  in  the  visibil¬ 
ities,  usually  represented  as  a  nontrivial  visibility  phase.  Recently 
developed  techniques  for  recovering  a  closely  related  quantity, 
“closure  phase”  (Baldwin  et  al.  1996;  Benson  et  al.  1997),  are 
now  being  applied  to  the  first  round  of  stellar  objects,  (e.g., 
Wittkowski  et  al.  2001). 

Although  originally  proposed  as  a  follow-up  on  the  oblateness 
observations,  Altair  was  observed  at  the  Navy  Prototype  Optical 
Interferometer  (NPOI;  Armstrong  et  al.  1998)  while  the  three- 
beam  combiner  was  in  operation,  allowing  measurement  of  closure 
phase  around  one  complete  triangle.  Examination  of  the  data  im¬ 
mediately  revealed  the  intermediate  phase  angles,  unambiguously 
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TABLE  1 

Roche  Model  Fits  for  Altair  and  (  Aql 


Altair  f  Aql 


Parameter  V\ 23  and  0c|a  Errors  Vf  and  tf>c\ b  /3  =  0.09  Vf  and  <pc\ b  Errors 

(1)  (2)  (3)  (4)  (5)  (6)  (7) 


OJ  =  n/Qc .  0.90  ±0.02  0.88  0.978  0.990  ±0.005 

6p  (mas) .  2.96  0.04  3.00  3.04  0.815  0.005 

Tp  (K) .  8740  140  8600  7980  11750 

/(deg) .  63.9  1.7  62.2  65.6  90  +0,-5 

P.A.  (deg) .  123.2  2.8  120.4  97.4  45  ±5 

X2/dof .  3.8  ...  9.5  13.4 

V  (obs:  0.77) .  0.765  ...  0.765  0.76 

/>’  /  lobs:  0.22) .  0.215  ...  0.22  0.26 

Veq  sin  i  ( km  s-1) .  245  ...  231  295 


a  Model  fit  to  the  triple  amplitude  and  closure  phase  data. 

b  Models  fit  to  the  closure  phase  data  and  the  squared  visibilities  of  the  three  baselines. 


signaling  the  presence  of  an  asymmetric  intensity  profile  (Ohishi 
et  al.  2003). 

Using  a  model  consisting  of  a  limb-darkened  disk  and  a  bright 
spot,  Ohishi  et  al.  (2004)  demonstrated  both  the  previously  dis¬ 
covered  oblateness  and  the  necessity  of  including  asymmetries 
in  the  intensity  distribution.  They  argued  that  the  probable  inter¬ 
pretation  was  that  of  rotational  flattening  and  gravity  darkening. 

In  the  meantime,  we  have  become  aware  of  some  limitations 
in  those  data  due  to  inadequate  corrections  for  “deadtimes”  in 
the  avalanche  photodiode  detectors  that  affect  the  high  signal 
levels  from  objects  as  bright  as  Altair.  We  therefore  reconsider 
a  subset  of  these  data  that  is  relatively  immune  to  the  detector 
problems,  using  a  full  implementation  of  the  von  Zeipel  theory 
(von  Zeipel  1924)  for  the  model  fitting  and  redoing  the  reduc¬ 
tions  in  a  way  that  dramatically  reduces  noise  in  the  bluest  chan¬ 
nels.  We  find  that  a  Roche  model  rotating  at  90%  of  the  breakup 
angular  velocity  and  inclined  ^64°  from  pole-on  fits  the  obser¬ 
vations  with  high  fidelity. 

We  show  that  the  parameter  that  sets  the  overall  temperature 
scaling  for  the  model,  the  effective  temperature  at  the  poles,  Tp, 
is  close  to  8700  K  for  this  model,  and  the  polar  surface  gravity 
is  correspondingly  fairly  high.  This  suggests  that  Altair  is  less 
evolved  than  one  might  naively  expect  from  its  spectral  type  and 
luminosity  classifications. 

In  this  model  the  equator  is  1 850  K  cooler  than  the  pole.  Given 
that  the  model  includes  both  polar  brightening  and  a  long  equatorial 
swath  of  low  intensity,  this  is  a  complex  intensity  distribution,  and 
the  agreement  with  the  observations  is  a  strong  endorsement  for 
the  simple  von  Zeipel  (1924)  theory. 

Below  we  describe  the  new  reductions,  give  a  brief  review  of 
Roche  theory,  and  then  present  the  fits.  We  note  that  the  exis¬ 
tence  of  large  amounts  of  surface  at  near-solar  temperatures  sug¬ 
gests  that  the  role  of  Altair  (and  probably  a  Cep  )  in  defining  the 
high-temperature  end  of  convection  on  the  main  sequence  may 
need  to  be  reconsidered.  We  also  note  the  recent  announcement 
that  Altair  is  a  low-amplitude  ( 6  Set)  pulsating  star  (Buzasi  et  al. 
2005),  which  may  give  hope  that  asteroseismology  will  be  able 
to  put  useful  limits  on  any  gradient  of  the  angular  velocity  in  the 
outer  envelope. 

2.  OBSERVATIONS 

Altair  was  observed  on  four  nights,  2001  May  25-27  and 
June  1 ,  with  the  NPOI.  These  are  the  same  observations  used  by 
Ohishi  et  al.  (2003, 2004);  we  refer  the  reader  to  those  papers  for 
a  journal  of  observations  and  a  description  of  the  observing  de¬ 


tails,  but  we  briefly  reprise  them  here.  We  have  focused  here  on 
the  data  set  obtained  2001  May  25.  This  is  by  far  the  largest  set  of 
data,  while  the  other  data  do  not  increase  the  range  of  hour  angles 
observed  in  the  first  night. 

The  observations  used  the  Astrometric  West  (AW),  Astrome¬ 
tric  East  (AE),  and  West  7  (W7)  stations,  forming  a  triangle  of 
interferometric  baselines  with  lengths  of  37.5  m  (AW-AE),  29.5  m 
(W7-AW),  and  64.4  m  (AE-W7).  The  back  end  combined  these 
three  input  beams  to  produce  three  output  beams,  with  one  base¬ 
line  on  each.  The  output  beams  were  dispersed  into  32  spectral 
channels  covering  22443 — 852  nm,  although  the  bluest  four  chan¬ 
nels  (22443— 460  nm)  of  the  W7-AW  output  were  not  functioning. 

The  Altair  observations  were  interleaved  with  observations 
of  a  calibrator,  £  Aql  (A0  V),  about  12°  away  on  the  sky.  We  ini¬ 
tially  estimated  its  diameter  to  be  around  0.85  mas,  with  which  it 
would  have  acted  as  a  quite  acceptable  calibrator.  However,  as 
noted  by  Ohishi  et  al.  (2004),  £  Aql  is  a  rapid  rotator  with  values 
of345  km  s”1  (Uesugi  &  Fukuda  1982)  and  317  km  s_1  (Royer 
et  al.  2002)  reported.  We  adopted  325  km  s  " 1 ,  raising  the  ques¬ 
tion  of  aspect  dependent  corrections  to  the  squared  visibilities 
and  phases.  This  possibility  was  discussed  by  Ohishi  et  al.  (2004), 
who  concluded  that  the  effects  were  not  important  at  the  level  of 
the  analysis  they  conducted,  but  cautioned  that  the  problem  needed 
to  be  reconsidered  if  a  detailed  analysis  of  Altair  was  attempted 
with  these  observations. 

We  have  found  a  number  of  occasions  in  2004  when  £  Aql 
was  observed  with  a  second  calibrator,  7  Lyr  (B9  111).  7  Lyr  is 
a  relatively  slow  rotator  (70  km  s  1 ;  Royer  et  al.  2002)  and  a  bit 
fainter  (V  ~  3.24),  which,  coupled  with  slightly  higher  temper¬ 
ature,  leads  to  the  expectation  of  a  symmetric,  nearly  unresolved 
calibrator  for  £  Aql.  This  turned  out  to  be  correct,  and  to  our  sur¬ 
prise,  we  found  that  we  were  able  to  deduce  a  meaningful  fit  of 
a  Roche  model  to  £  Aql  using  the  phase  and  visibility  amplitude 
data.  Since  such  results  are  rare  (our  report  on  Altair  here  is 
the  first),  we  have  decided  to  present  a  detailed  discussion  of  the 
case  in  a  separate  communication  (D.  M.  Peterson  et  al.  2006,  in 
preparation). 

We  summarize  in  columns  (6)  and  (7)  of  Table  1  a  prelimi¬ 
nary  set  of  the  relevant  Roche  parameters  for  £  Aql.  Given  the 
small  angular  diameter,  the  quoted  errors  produce  uncertainties 
in  the  calibration  of  Altair  that  are  undetectable  compared  to 
other  error  sources.  One  property  of  note  from  these  parameters 
is  the  near  equator-on  orientation  of  the  spheroid.  At  most,  we 
would  have  expected  a  few  degrees  uncertainty  in  the  triple 
phases  induced  by  the  range  of  possible  inclinations.  However, 
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the  observations  calibrated  by  7  Lyr  show  that  these  phases  are 
truly  small,  indistinguishable  from  zero  at  the  ±  1  °  level.  The  star 
is  apparently  seen  at  nearly  i  ~  90°. 

Ohishi  et  al.  (2004)  used  contemporaneous  observations  of 
Vega  as  a  check  star — one  that  should  show  a  circular  outline.  Al¬ 
though  we  agree  those  observations  do  seem  to  imply  that  Vega 
is  circular,  we  are  not  inclined  to  place  much  weight  on  the  re¬ 
sult.  Vega  is  twice  as  bright  as  Altair  in  this  wavelength  range. 
Detector  nonlinearities,  already  serious  in  the  Altair  observations 
(as  we  describe  in  §  3),  are  overwhelming  here.  We  simply  do  not 
know  how  to  interpret  the  observed  squared  visibilities. 

Finally,  we  note  that  since  the  NPOI  records  phases  and  vis¬ 
ibilities  over  a  range  of  wavelengths,  there  is  little  possibility 
of  contamination  by  unknown  companions,  at  least  those  within 
the  0ff5  field  of  view  of  the  siderostats.  If  such  existed,  there 
would  be  clear,  strongly  modulated  phases  and  visibilities.  (  Aql 
is  a  known  member  of  a  wide  multiple  system,  ADS  12026,  with 
companions  no  closer  than  5"  and  magnitude  differences  no  less 
than  8. 

3.  DATA  REDUCTION 

3.1.  Incoherent  Integration 

The  NPOI  observes  interference  fringes  by  modulating  the 
optical  path  on  the  delay  line  for  each  array  element,  using  a 
triangle-wave  pattern  at  a  frequency  of  500  Hz.  The  resulting 
modulation  of  the  intensity  is  detected  in  eight  bins  evenly  spaced 
over  one  fringe  in  each  channel  by  avalanche  photodiodes.  The 
phase  of  the  intensity  modulation  changes  on  timescales  of  milli¬ 
seconds  since  the  fringe  tracker  employed  by  NPOI  tracks  the  en¬ 
velope  of  the  (bandwidth-limited)  fringe  packet  rather  than  the 
fringe  phase. 

The  data  in  the  delay  bins  were  processed  to  produce  the  com¬ 
plex  visibility  and  squared  visibility  modulus  V2  for  each  base¬ 
line  at  each  wavelength.  From  these,  the  triple  product  K123  exp  i(pc\ 
can  be  calculated,  where  the  triple  amplitude  F123  =  |  V\  1 1 V2 1 1 V3 1 
is  the  product  of  the  amplitudes  of  the  complex  visibilities  of  the 
individual  baselines  and  the  closure  phase  <j>c \  =  <j>\  +  (j>i  +  </>3  is 
the  sum  of  the  individual  phases.  Although  the  baseline  phases 
themselves  are  affected  by  atmospheric  turbulence,  those  effects 
cancel  in  the  sum  of  three  phases  around  a  closed  triangle,  so  the 
closure  phase  preserves  information  about  the  source  structure. 
These  data  products  are  produced  for  each  2  ms  cycle  of  delay 
modulation.  In  the  standard  incoherent  integration  as  described 
by  Hummel  et  al.  (1998),  the  squared  visibilities  and  complex 
triple  products  are  summed  to  provide  average  values  in  1  s 
intervals. 

3.2.  Coherent  Integration 

We  employed  a  new  algorithm  for  the  coherent  integration  of 
the  complex  visibilities  of  the  NPOI  first  presented  by  Hummel 
et  al.  (2003).  Compared  to  the  incoherent  integration  of  the 
squared  visibilities,  coherent  integration  achieves  a  higher  signal- 
to-noise  ratio  of  the  averages  due  to  the  larger  number  of  photons 
detected  in  a  coherent  sample  of  the  fringe.  We  have  exploited 
this  fact  to  recover  meaningful  results  from  all  NPOI  spectrom¬ 
eter  channels,  while  the  channels  on  the  blue  side  of  about  560  nm 
had  usually  been  discarded  in  incoherent  reductions  due  to  the 
insufficient  number  of  photons  detected  during  a  2  ms  instrumen¬ 
tal  integration  time.  For  the  coherent  integration  time  we  selected 
200  ms,  and  the  resulting  complex  visibilities  were  both  com¬ 
bined  to  form  complex  triple  products  and  transformed  individ¬ 
ually  into  squared  amplitudes  of  the  modulus.  Every  1 0  samples  of 
these  quantities  were  then  averaged  (averaging  real  and  imaginary 


part  of  the  complex  triple  products  separately)  for  a  total  integra¬ 
tion  time  per  data  point  of  2  s. 

The  alignment  of  the  raw  visibility  phasors,  necessary  be¬ 
fore  integration  in  order  to  avoid  detrimental  coherence  losses, 
was  performed  as  follows.  Two  steps  are  necessary  to  rotate  the 
phasors  onto  a  common  fringe  in  order  to  enable  a  phase-tracking 
algorithm. 

First,  we  computed  the  average  power  spectra  of  the  chan¬ 
neled  visibility  as  a  function  of  delay  for  10  ms  intervals.  Their 
maxima,  corresponding  to  the  group  delays  of  the  fringe  packets, 
are  not  zero  but  have  a  typical  mis  on  the  order  of  1  gm  as  the 
NPOI  group  delay  fringe  tracker  tries  to  center  the  fringe  but 
does  not  lock  onto  its  phase. 

Second,  from  the  deviation  in  position  of  the  fringe  from  the 
estimated  geometrical  value,  which  is  nonzero  and  has  a  typical 
mis  on  the  order  of  10  gm  due  to  atmospheric  refractive  index 
fluctuations,  we  estimated  the  differential  amount  of  air  and  thus 
the  phase  shift  between  the  peak  of  the  envelope  of  the  fringe 
packet  and  the  nearest  fringe  peak.  In  other  words,  this  phase  is 
the  phase  of  the  complex  Fourier  transform  of  the  visibility  as  a 
function  of  wavenumber.  The  modulus  of  this  transform  peaks  at 
the  value  of  the  group  delay.  The  phase  of  the  transform  at  this 
delay  is  called  the  group  delay  phase. 

We  converted  the  group  delay  phase  to  a  delay  using  the  mean 
wavelength  of  the  white-light  fringe  and  added  it  to  the  group 
delay.  Rotation  of  the  visibility  phasors  of  different  channels  by 
an  angle  corresponding  to  the  ratio  of  this  delay  value  and  the 
wavelength  of  the  spectrometer  channel  will  align  them  on  the 
same  fringe.  At  this  point,  the  algorithm  implements  a  photon- 
noise  limited  off-line  fringe  phase  tracker  enabling  the  use  of 
much  longer  coherent  integration  times. 

3.3.  Baseline  Bootstrapping 

We  used  an  important  modification  of  the  above  procedure  by 
applying  the  baseline  bootstrapping  method,  a  design  feature  of 
the  NPOI  interferometer  (Armstrong  et  al.  1998).  It  exploits  the 
fact  that  the  sum  of  the  fringe  delays  along  a  closed  loop  of  base¬ 
lines  is  zero  (if  the  same  fringe  is  identified  on  each  baseline). 
Therefore,  if  a  long  baseline  in  a  multitelescope  array  sees  a  low 
contrast  fringe  due  to,  for  example,  object  extension,  and  this 
baseline  involves  two  telescopes  that  are  at  the  same  time  in¬ 
volved  with  other  telescopes  of  the  array  on  much  shorter  base¬ 
lines  seeing  much  higher  fringe  contrast,  the  fringe  delay  of 
the  long  baseline  can  be  computed  from  the  fringe  delays  on  the 
shorter  baselines,  which  are  “bootstrapping”  the  long  one.  In  the 
simple  case  of  the  observations  described  on  Altair,  the  fringe 
delay  on  the  long  64  m  W7-AE  baseline  is  just  the  difference 
between  the  fringe  delays  on  the  shorter  AE-AW  and  W7-AW 
baselines. 

3.4.  Averaging  and  Editing 

The  2  s  data  points  are  edited  for  outliers  as  described  by 
Hummel  et  al.  (1998).  The  final  averaging  is  done  over  the  full 
length  of  a  pointing  (  called  a  scan  at  NPOI  ),  which  typically  lasts 
90  s.  The  computation  of  the  formal  errors  also  follows  Hummel 
et  al.  (1998),  except  that  we  have  implemented  a  different  ap¬ 
proach  for  the  complex  triple  products  based  on  a  suggestion  by 
D.  Buscher  (2003,  private  communication).  Under  simple  as¬ 
sumptions,  the  error  of  a  complex  triple  product  is  described  by 
an  error  ellipse  that  has  one  axis  aligned  with  the  triple  product 
phasor.  Assuming  this,  we  compute  the  error  of  the  triple  product 
as  the  error  of  the  imaginary  and  real  parts  of  the  mean  after  ap¬ 
plying  a  rotation  of  all  phasors  by  the  mean  triple  product  phasor. 
The  errors  of  amplitude  and  phase  of  the  triple  product  are  then 
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equal  to  the  error  of  the  real  part  and  the  error  of  the  imaginary 
part  divided  by  the  respective  amplitudes. 

3.5.  Detector  Nonlinearity’ 

A  source  of  systematic  error  comes  from  deadtime  in  the  pulse 
counting  electronics  controlling  the  avalanche  photodiode  de¬ 
tectors.  These  systems  saturate  at  about  1  MHz  and  display  sig¬ 
nificant  nonlinearity  in  the  apparent  count  rates  as  they  near 
this  limit.  The  nominal  design  of  the  detector  systems  included  a 
r  =  200  ns  deadtime,  but  we  have  subsequently  found  that  not 
only  do  those  time  constants  vary  significantly,  channel  to  chan¬ 
nel,  but  also  that  in  a  given  channel  they  depend  on  the  mean 
signal  level  because  of  the  effects  of  heating.  We  believe  that  it 
will  be  possible  to  model  and  remove  these  effects,  but  some  ef¬ 
fort  is  involved,  which  we  will  report  on  in  the  future.  Unfor¬ 
tunately,  these  problems,  which  do  not  affect  the  fainter  objects 
usually  observed  by  the  instrument,  were  not  recognized  at  the 
time  the  Altair  observations  were  made. 

However,  we  believe  that  through  a  rather  unique  set  of  cir¬ 
cumstances  the  phase  and  some  of  the  amplitude  measurements 
acquired  during  the  2001  observations  are  to  first  order  free  of 
the  effects  of  these  nonlinearities.  One  reason  was  that  during 
these  observations  only  three  stations  were  in  use,  and  the  three 
spectrographs  recorded  single  baseline  data. 

The  other  reason  was  that  the  amplitude  and  phase  measure¬ 
ments  from  each  channel  were  accomplished  with  a  simple  dis¬ 
crete  Fourier  transform  (DFT).  By  dithering  the  optical  delay  at 
a  frequency  of  u>,  the  signal  was  modulated  according  to 

I(t)  =/0[l  +  V  cos(ojt  +  </>)],  (1) 

where  V  is  the  (instrumental)  amplitude  of  the  visibility  and  (j> 
the  instantaneous  phase.  The  detector  system  responded  to  the 
modulated  signal  according  to 

N(0  =  l  (f^Qj{t)  ~  *o(0  -  rNl{t)  +...,  (2) 

where  Q  is  the  quantum  efficiency,  N(t)  is  the  apparent  pulse  rate, 
and  No(t)  =  QI(t )  is  the  true  photon  detection  rate.  In  the  line¬ 
arized  form  we  assumed  tNq  1 ,  which  for  a  300  kHz  count  rate, 
typical  of  the  wider  red  channels  on  the  higher  visibility  base¬ 
lines,  is  adequate  to  1%  or  better.  By  substituting  equation  (1) 
into  equation  (2),  clearing  the  quadratic  cosine  using  the  half¬ 
angle  formula,  and  collecting  terms,  this  becomes 

N(t)  ~  No  [  1  -  tN0  ( 1  +  F2/2)  ]  +  N0  V(  1  -  2rN0)  cos  (cct  +  </>) 
tNSV2 

- —  cos  (2c et  +  2(f>),  (3) 

where  No  =  QIq.  A  DFT  at  the  dither  frequency  now  extracts 
an  amplitude  different  than  the  nominal  No  V  (and  after  division 
by  the  nominal  mean  signal  produces  an  estimate  that  can  differ 
significantly  from  the  true  visibility  amplitude).  However,  the 
phase  comes  through  the  process  unaffected,  as  do  the  frequen¬ 
cies  of  the  minima  in  the  visibility  amplitudes.  By  floating  the 
overall  amplitudes  in  the  data  reduction  (see  below),  we  retain 
the  important  spatial  scale  information  contained  in  the  minima. 
But  the  most  important  conclusion  is  that  the  phases  may  be  as¬ 
sumed  to  be  essentially  free  of  detector  induced  biases. 

3.6.  Visibility >  Calibration 

The  degradation  of  the  measured  visibilities  due  to  atmo¬ 
spheric  and  instrumental  effects  is  measured,  as  with  all  inter¬ 


ferometers,  by  observing  calibrator  stars  with  diameters  as  small 
as  possible  to  reduce  uncertainties  in  the  visibility  estimates  for 
them.  As  shown  by  Hummel  et  al.  (1998),  the  NPOI  visibility 
amplitudes  sometimes  show  a  negative  correlation  with  the  rms 
of  the  delay  line  motion  that  is  related  to  the  seeing.  But  at  other 
times,  instrumental  effects  that  correlate  with  time  or  other  sys¬ 
tematic  effects  that  correlate  with  position  on  the  sky  (e.g.,  hour 
angle)  can  dominate  the  visibility  variations.  Therefore,  formal 
photon  noise-based  visibility  errors  usually  require  the  addition 
(in  quadrature)  of  a  calibration  error  that  is  derived  from  the  re¬ 
sidual  visibility  variations  of  the  calibrator  after  calibration.  For 
the  amplitude  calibrations,  we  smoothed  the  calibrator  visibili¬ 
ties  with  a  20  minute  Gaussian  kernel  in  hour  angle  and  obtained 
calibration  errors  ranging  from  about  4%  at  the  red  end  to  1 5%  at 
the  blue  end  of  the  spectrometers.  For  the  closure  phase,  we  used 
the  same  smoothing  technique  but  applied  to  the  calibrator  phases 
as  a  function  of  time. 

During  model  fitting,  however,  this  standard  procedure  yielded 
inconsistent  fits.  One  can  already  see  this  from  the  results  of 
Ohishi  et  al.  (2004,  their  Fig.  6),  where  amplitudes  can  be  sys¬ 
tematically  high  or  low  with  the  important  characteristic  that  the 
deviation  is  very  consistently  independent  of  wavelength.  The 
effect  is  exacerbated  by  the  fact  that  due  to  the  brightness  of 
Altair,  the  formal  amplitude  errors  are  quite  small.  The  reason  for 
the  scan-to-scan  variations  is  most  likely  the  same  as  for  the 
residual  variations  of  the  calibrator  after  calibration,  except  that 
there  is  no  perfect  correlation  due  to  target  and  calibrator  not 
being  at  the  same  location  in  the  sky.  (Past  experience  has  shown 
that  visibilities  do  correlate  quite  well  if  the  calibrator  is  very 
near  the  target.)  Therefore,  we  allowed  the  calibration  for  each 
scan  and  baseline  to  float  by  applying  “achromatic”  calibration 
factors  to  improve  the  fit  between  data  and  model.  We  discuss  the 
implications  for  the  model  fitting  in  §  4.  Finally,  we  note  that  an 
error  in  the  calibrator  star  diameter  produces  chromatic  errors 
across  the  NPOI  spectrometer,  as  do  uncompensated  deadtime 
corrections,  neither  of  which  can  be  removed  by  the  calibration 
factors. 

4.  MODELING 
4.1.  Roche  Spheroids 

The  theory  for  the  equilibrium  shapes  and  surface  properties 
of  rotating  stars  was  first  presented  80  years  ago  (von  Zeipel 
1924)  assuming  solid  body  rotation  and  a  point  source  gravita¬ 
tional  potential.  This  model  has  proved  quite  successful  in  de¬ 
scribing  the  figures  of  stars  in  close  binary  systems,  where  tidal 
effects  to  first  order  produce  the  same  distortions  as  rotation 
(Collins  1989). 

We  use  that  model  here  but  note  its  limitations.  First,  there  is 
no  a  priori  reason  that  stars  should  rotate  as  solid  bodies,  and  the 
surface  layers  of  the  Sun  have  long  been  known  to  rotate  dif¬ 
ferentially.  However,  among  the  early  results  of  helioseismology 
was  the  discovery  that  the  transition  from  the  outer  convection 
zone  to  the  inner  radiative  layers  coincided  with  an  abrupt  tran¬ 
sition  to  solid  body  rotation  (Spiegel  &  Zahn  1992).  Since  early- 
type  stars  have  radiative  envelopes  and  relatively  small  (in  radius) 
convective  cores,  one  might  expect  solid-body  rotation  to  be  a 
good  approximation  for  the  external  layers  of  early-type  stars. 

There  has  been  some  observational  support  for  this  expecta¬ 
tion.  Reiners  &  Royer  (2004)  have  analyzed  the  rotational  profiles 
of  a  large  number  of  A  stars  looking  for  evidence  of  differen¬ 
tial  rotation  following  a  solar-type  latitudinal  dependence.  In  the 
78  stars  for  which  the  determination  could  be  made,  they  found 
four  objects  for  which  peculiarities  were  seen  that  might  be  from 
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differential  rotation  (or  other  causes).  However,  95%  of  the  line 
shapes  were  fully  consistent  with  solid-body  rotation. 

In  addition,  there  is  the  long-known  consistency  between  the 
largest  rotational  velocities  measured  in  the  early-type  stars  and 
the  predicted  maximum  rotation  velocities  associated  with  “equa¬ 
torial  breakup”  (Fremat  et  al.  2005).  In  recent  decades  (see,  e.g., 
Tassoul  1978)  it  has  been  demonstrated  that  rotation  laws  other 
than  rigid  rotation  do  not  generally  impose  maximum  rotation 
velocities. 

The  second  limitation,  one  that  we  will  spend  some  time  on, 
involves  the  exponent  in  the  TeS  —  geS  relation  (geff,  the  effec¬ 
tive  gravity,  includes  centrifugal  terms).  In  the  original  work,  von 
Zeipel  (1924)  considered  the  case  of  a  fully  radiative  envelope, 
deriving  the  well-known  “gravity  darkening”  relation  Teg  <x 
g°g5.  Lucy  (1967)  reconsidered  the  problem  in  fully  convective 
stars,  deriving  a  much  reduced  gravity  dependence,  Tes  ex  g^8. 
Other  approximations  lead  to  other  exponents  (see  Reiners  2003 
for  references).  In  our  nominal  calculations  we  adopt  the  original 
von  Zeipel  (1924)  prescription.  As  we  show  below,  the  Altair  ob¬ 
servations  bear  significantly  on  this  issue. 

Even  in  the  limit  of  rigid  rotation,  the  von  Zeipel  (1924) 
theory  is  only  first  order  in  rotation  rate.  Distortions  of  the  in¬ 
terior  figure,  allowing  some  gravitational  quadrapole  contribu¬ 
tion,  can  be  expected  as  rotation  rates  approach  breakup.  That  is 
even  more  likely  if  there  are  significant  deviations  from  solid- 
body  rotation,  even  if  confined  to  the  inner  convective  regions 
(e.g.,  Tassoul  1978). 

Probably  even  more  relevant  are  the  effects  of  radiation  pres¬ 
sure,  which  is  treated  simplistically  in  the  theory,  and  stellar  winds. 
Significant  envelope  extension  due  to  radiation  can  be  expected  in 
the  low  effective  surface  gravity  regions  of  rapidly  rotating  stars. 
And  we  would  certainly  expect  a  dramatic  increase  in  mass  loss  at 
the  equator,  both  due  to  enhanced  convection  as  gas  temperatures 
decrease  into  the  solar  and  subsolar  regime.  We  ignore  all  these 
effects  here,  except  to  acknowledge  the  limitations  inherent  in  this 
theory. 

4.2.  Roche  Models 

We  have  constructed  a  suite  of  programs  to  evaluate  the  run 
of  specific  intensity  across  the  surface  of  a  Roche  spheroid.  The 
definitions  of  the  various  angles  are  from  Collins  (1963).  Other¬ 
wise,  we  follow  the  prescription  for  the  surface  figure  and  the 
notation  given  by  Hardorp  &  Strittmatter  (1968)  with  one  ex¬ 
ception:  following  the  discussion  by  Hardorp  &  Scholz  (1971) 
we  take  the  polar  radius,  Rp,  as  a  fixed  parameter.  Specifically, 
we  do  not  allow  it  to  be  a  function  of  the  fractional  rotation. 

The  modeling  requires  that  we  specify  six  quantities:  the  ratio 
of  the  angular  velocity  to  that  of  breakup,  lu  =  Cl/Clg,  the  in¬ 
clination  (or  tilt)  of  the  rotational  axis,  i,  defined  such  that  i  =  0 
is  pole-on,  the  position  angle,  P.A.,  of  the  pole  on  the  sky  (mea¬ 
sured  north  through  east),  the  angular  diameter  of  the  polar  axis, 
dp,  the  effective  temperature  at  the  pole,  Tp ,  and  the  surface  grav¬ 
ity,  or  more  commonly  the  logarithm  of  the  surface  gravity  (cgs), 
at  the  pole,  log  gp.  From  the  relations  in  the  cited  references  it  is 
then  possible  to  calculate  the  radius  of  the  star  for  a  given  stellar 
latitude,  9,  and  hence  the  surface  gravity  geff  and  effective  tem¬ 
perature  [ Ta(6 )  =  Tp(geS/gp)\  at  that  latitude  (see  Hardorp  & 
Strittmatter  1968  for  the  definition  of  geff). 

Note  that  for  stars  with  accurate  parallaxes  like  Altair,  spec¬ 
ifying  the  polar  surface  gravity  and  angular  diameter  is  equiv¬ 
alent  to  specifying  the  mass  and  linear  (polar)  radius.  These  in 
turn  fix  the  breakup  angular  velocity,  VL2b  =  (8/27 )GM/R2,  along 
with  the  equatorial  and  projected  velocities  (when  the  inclination 
is  specified).  Finally,  it  is  useful  to  recall  the  relation  between  the 


polar  radius  and  equatorial  radius  at  breakup  angular  velocity: 
Re  b  =  3Rp/2  (  Hardorp  &  Strittmatter  1968).  According  to  this 
first-order  theory,  the  maximum  rotational  flattening  is  |. 

The  model  definition  is  completed  by  specifying  for  each 
wavelength  and  surface  point  the  specific  intensity  at  the  angle 
of  the  line  of  sight  from  the  local  normal.  As  noted  by  the  early 
authors,  plane-parallel  model  atmospheres  are  entirely  adequate 
in  the  context  of  Roche  models  for  stars  on  and  near  the  main 
sequence.  The  only  exception  is  that  these  models  develop  a 
cusp  at  the  equator  at  critical  rotation  velocity.  However,  the  cusp 
does  not  appear  until  fractional  rotation  velocities  of  lo  =  0.99  or 
larger,  and,  as  we  have  indicated,  the  breakdown  of  the  plane- 
parallel  approximation  is  only  one  of  several  problems  with  the 
model  in  this  limit.  A  number  of  auxiliary  quantities  are  also  cal¬ 
culated,  both  as  integrated  over  the  spheroid  and  as  seen  from  a 
particular  direction,  in  addition  to  the  predicted  complex  visibil¬ 
ities.  An  example  of  the  former  is  the  integrated  luminosity,  cal¬ 
culated  as  described  in  Collins  (1963).  Examples  of  the  latter 
include  the  magnitudes  in  various  bandpasses. 

In  practice,  we  need  to  solve  the  inverse  problem:  given  a 
point  on  the  sky  (a,  6),  determine  whether  the  point  is  on  the 
stellar  disk  and  if  so  what  the  corresponding  latitude  and  radius 
are.  We  have  solved  this  problem  explicitly  using  simple  itera¬ 
tion  and  first-  and  second-order  versions  of  Newton-Raphson 
iteration.  The  routines,  written  in  C,  are  quite  flexible,  reason¬ 
ably  fast,  and  freely  available  (from  the  first  author). 

The  properties  of  Roche  models  for  isolated  rotating  early- 
type  stars  have  recently  been  reviewed  by  Domiciano  de  Souza 
et  al.  (2002),  who  use  a  slightly  different  but  completely  equiv¬ 
alent  parameterization.  Those  authors  summarize  some  of  the 
results  from  their  models  for  massive  stars,  which  have  provided 
a  useful  check  on  our  own  routines. 

4.3.  Implementation  of  Roche  Code 

The  Roche  code  consists  of  a  library  of  functions  written  in  C, 
with  a  main  function  enabling  its  use  as  standalone  software  and 
a  wrapper  enabling  it  to  be  called  from  within  the  NPOI  standard 
data  reduction  software  OYSTER.  The  Roche  spheroid  parame¬ 
ters  are  part  of  the  standard  hierarchical  model  format  of  O  Y  STER 
and  are  passed  to  the  Roche  code  along  with  pointers  to  the  ex¬ 
tensive  tables  of  linear,  logarithmic,  and  square-root  law  mono¬ 
chromatic  limb-darkening  coefficients  for  a  grid  of  Kurucz  model 
atmospheres  as  published  by  Van  Hamme  ( 1 993).  The  Roche  code, 
with  the  additional  input  of  the  ( u ,  v)  coordinates,  computes  the 
visibilities  for  a  grid  of  wavelengths  supplied  by  OYSTER,  which 
are  subsequently  integrated  over  the  NPOI  bandpasses. 

The  fitting  of  the  model  parameters  (except  for  the  gravities) 
utilizes  the  Marquardt-Levenberg  algorithm  (Press  et  al.  1992) 
implemented  in  OYSTER,  with  the  derivatives  computed  nu¬ 
merically.  In  addition  to  the  visibility  and  phase  measurements, 
the  reductions  were  also  constrained  to  reproduce  the  observed  V 
magnitude.  This  provided  particularly  strong  constraints  on  the 
polar  effective  temperature. 

5.  DISCUSSION 
5.1.  Model  Fitting 

The  analysis  proceeded  with  few  complications.  In  particular, 
a  close  examination  of  the  y2  surface  indicated  no  unusual  mor¬ 
phology,  and  indeed  the  iterations  converged  to  the  same  final 
solution  independent  of  our  starting  guess,  whether  from  larger 
or  smaller  values  of  the  parameters.  Our  final  solution,  given  in 
column  (2)  of  Table  1,  is  based  on  the  triple  phases  and  triple 
amplitudes  only,  the  latter  with  an  overall  floating  multiplier  for 
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Fig.  1 . — Triple  amplitudes  as  fit  by  best  Roche  model.  Error  bars  of  ±cr  are  shown.  The  analytic  fits  ( solid  lines)  include  a  constant  multiplicative  renormalization. 


each  scan,  as  described  in  §  3.6.  The  reduced  x2  for  this  solu¬ 
tion  is  still  a  bit  large,  which  we  attribute  to  some  residual  non¬ 
gray  problems  of  uncertain  origin  in  the  amplitudes,  as  shown  in 
Figure  1 .  The  fit  to  the  closure  phases  (Fig.  2),  on  the  other  hand, 
is  remarkable,  showing  no  trends  with  wavelength  or  hour  angle. 

The  for  our  best  model  is  well  in  excess  of  unity.  This  is 
the  result  of  both  the  remaining  residuals  in  the  triple  amplitudes 
evident  in  Figure  1  and  the  very  small  formal  errors  in  both 
the  amplitudes  and  phases  owing  to  the  high  signal  levels.  If  the 
large  x2  had  resulted  from  just  the  latter,  we  could  have  taken  the 
usual  expedient  and  scaled  the  errors  by  a  common  multiplier, 
which  would  have  normalized  away  the  excess  x2-  Estimates  of 
the  uncertainties  in  the  resulting  parameters  would  then  have 
been  obtained  through  the  usual  process  of  varying  each  param¬ 
eter  until  x2  had  increased  by  unity,  thereby  mapping  out  formal 
errors  and  the  correlation  matrix. 

However,  this  gives  very  small  estimated  errors  that  we  feel  do 
not  properly  reflect  the  influence  of  the  (small)  remaining  biases 
in  the  triple  amplitudes.  We  have  therefore  decided  to  take  a  con¬ 
servative  approach  to  our  error  estimates.  We  have  run  a  separate 
reduction  using  the  squared  visibilities  for  the  individual  base¬ 
lines  in  place  of  the  triple  amplitudes.  This  solution,  which  has  a 
substantially  larger  reduced  x2  reflecting  the  larger  residuals  in 
the  amplitudes  of  the  individual  baselines,  is  summarized  in  col¬ 


umn  (4)  of  Table  1 .  Our  adopted  errors,  shown  in  column  (3),  are 
the  difference  between  these  two  solutions. 

Figure  3  shows  how  Altair  appears  projected  on  the  sky.  The 
intensity  distribution  at  500  nm,  as  would  be  seen,  for  example, 
by  an  interferometer,  is  color-encoded:  blue  for  high  intensity, 
red  for  low  intensity.  Except  for  limb-darkening,  this  is  also  a 
temperature  encoding.  The  range  in  intensities,  a  factor  of  18,  is 
about  a  factor  of  2.5  more  than  would  be  expected  due  to  limb- 
darkening  alone  in  a  nonrotating  star  of  this  spectral  type. 

Included  in  Table  1  is  the  integrated  B  —  V  as  calculated  for 
the  models.  Although  we  force  all  models  to  reproduce  the  ob¬ 
served  V  magnitudes,  as  described  above,  this  constraint  does 
not  automatically  mean  B  —  V  will  be  reproduced.  The  fact  that 
the  color  does  agree  with  the  measurements  is  a  significant  con¬ 
sistency  check  on  the  models,  particularly  given  that  the  range 
of  temperatures  across  the  surface  could  produce  a  wide  range 
in  B  —  V .  A  number  of  ancillary  parameters  derived  from  the 
adopted  fit  are  given  in  Table  2  and  are  discussed  further  in  §  6. 

5.2.  Comparison  with  Previous  Results 

As  noted,  Altair  has  been  the  subject  of  a  number  of  inter¬ 
ferometric  measurements  over  the  years,  most  notably  by  the 
PTI  array  (van  Belle  et  al.  2001).  The  diameter  determinations 
are  not  simply  intercomparable,  even  if  we  consider  only  the 
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25-May-01  Triple  Phases 


Fig.  2. — Triple  phases  as  fit  by  the  best  Roche  model.  No  renormalization  has  been  applied.  Small  circles  are  used  for  the  observations  since  on  the  red  side  of  the 
plots  the  error  bars  tend  to  be  smaller  than  (and  fall  within)  the  width  of  the  line  showing  the  analytic  fit. 


projected  major  and  minor  axes  of  our  model,  since  the  PT1  ob¬ 
servations  were  fit  to  the  physical  dimensions  of  a  Roche  model, 
which  included  limb  darkening  but  not  gravity  darkening.  Still, 
our  determination  of  an  equatorial  radius  of  1.988  Rq  seems  in 
reasonable  accord  with  their  quoted  1.88  R@. 

More  problematic  are  the  reported  position  angles  of  the  ro¬ 
tational  pole.  First,  we  note  that  in  a  preliminary  report  of  this 
work  (Peterson  et  al.  2004)  the  pole  is  off  by  180°  due  to  a  sign 
error.  More  complicated  is  the  disagreement  between  the  posi¬ 
tion  angles  from  the  PTI  measurements,  —25°  ±  9°  compared  to 
our  123?2  ±  2?8.  Communication  with  G.  van  Belle  (2003,  pri¬ 
vate  communication)  indicated  that  the  (u,  v)  coordinates  of  each 
baseline  were  inadvertently  exchanged  in  their  analysis.  Fitting  a 
simple  unifomi  ellipse  model  to  the  PTI  data,  corrected  for  the 
component  swap,  yields  P.A.  =  \22°2,  in  excellent  agreement 
with  our  result. 

Recently,  Reiners  &  Royer  (2004)  reported  the  determination 
of  Altair’s  equatorial  rotational  velocity,  j;eq  <  245  km  s  *.  They 
analyzed  the  star’s  rotational  broadening  profile  to  determine  the 
first  two  zeros  of  its  Fourier  transform,  the  ratio  of  which  has 
been  shown  by  Reiners  (2003)  to  depend  on  the  equatorial  veloc¬ 
ity  rather  than  the  usual  v  sin  i.  This  is  a  new  approach  to  measur¬ 
ing  total  velocities  in  stars,  and  it  is  difficult  to  know  how  much 


weight  it  should  be  given.  One  notable  aspect  of  that  analysis 
was  the  adoption  of  an  exponent  for  the  gravity-darkening  law 
(/ 3  ~  0.09)  that  was  about  |  that  of  the  von  Zeipel  (1924)  value. 
We  discuss  this  aspect  of  calculating  rotationally  distorted  stars 
next. 

5.3.  Gravity  Darkening 

While  we  have  several  lines  of  reasoning  (described  in  §  4.1) 
that  lead  us  to  believe  that  solid-body  rotation  is  valid  for  these 
models,  the  situation  is  not  so  clear  with  regard  to  the  exponent 
on  the  gravity-darkening  law,  TeS  cx  gfff.  The  classical  work  of 
von  Zeipel  (1924)  would  seem  to  apply  to  an  A  star,  even  one 
somewhat  evolved,  yielding  [3  =  0.25.  However,  convection  does 
occur  in  late  A  stars,  particularly  in  the  photosphere  where  it 
competes  with  radiation  in  carrying  the  flux.  It  is  then  impor¬ 
tant  to  note  that  Lucy  (1967)  has  shown  that  for  small  distor¬ 
tions  the  appropriate  coefficient  in  fully  convective  envelopes  is 
closer  to  /?  ~  0.08. 

This  leaves  matters  in  a  somewhat  uncertain  state.  On  the  one 
hand,  there  are  stars  in  the  transition  region  between  having  fully 
convective  and  fully  radiative  envelopes,  Altair  arguably  one  of 
them,  and  there  is  no  obvious  guidance  in  choosing  an  appro¬ 
priate  value  for  this  parameter.  On  the  other  hand,  even  where  the 
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Altair  i=63.9 


RA  (mas) 

Fig.  3. — False-color  rendering  of  Altair’s  visible  surface.  Intensity  at  500  nm  increases  from  red  to  blue.  Except  for  the  effects  of  limb  darkening,  this  is  also  a 
map  of  temperature,  which  varies  from  8740  K  at  the  pole  to  6890  K  at  the  equator. 


envelopes  are  unambigiously  in  either  one  of  those  states  or  the 
other,  the  classical  results  apply  under  rather  different  circum¬ 
stances,  which  we  discuss  next. 

Both  the  von  Zeipel  (1924)  and  Lucy  (1967)  results  treat  rota¬ 
tion  as  a  perturbation.  However,  in  the  radiative  case  in  which 
uniform  rotation  is  adequate  (and  issues  like  mass  loss,  etc.  can 
be  ignored),  the  quantity  treated  as  a  perturbation  is  the  size  of 
the  quadrapole  moment  of  the  gravitational  field.  Since  stars  are 
centrally  condensed,  even  for  velocities  approaching  critical  the 
distortions  in  the  core  are  modest,  one  can  expect  that  the  analy¬ 
sis  given  by  von  Zeipel  (1924)  will  be  reasonably  accurate.  This 
has  been  found  to  be  true  in  practice  (Sackmann  1970). 

For  the  convective  case,  the  gravity  darkening  exponent  is  ob¬ 
tained  by  analyzing  the  adiabats  found  in  the  envelopes  of  rep¬ 
resentative  stars.  Lucy  (1967)  quite  explicitly  points  out  that  the 
derivation  is  valid  only  for  small  changes  in  the  effective  gravity. 


Of  course,  the  effective  gravity  changes  by  orders  of  magnitude 
as  rotation  approaches  critical,  and  it  is  not  clear  whether  the 
exponent  derived  by  Lucy  can  be  used  to  describe  gravity  dark¬ 
ening  for  anything  but  the  most  modest  rotation.  Again,  this  is  in 
contrast  to  the  modest  contribution  of  an  induced  gravitational 
quadrapole,  even  for  stars  rotating  at  breakup. 

Even  so,  in  a  recent  series  of  papers,  Claret  (2004  and  refer¬ 
ences  therein)  has  attempted  to  deal  with  the  issue  of  a  smooth 
interpolation  between  these  two  extreme  cases.  He  has  noted  that 
as  stars  evolve  off  the  main  sequence  and  toward  the  red  giant 
branch,  their  interior  structures  trace  out  approximately  straight 
line  loci  in  a  (log  T^,  log  g)  diagram.  On  the  main  sequence  for 
massive  (mostly  radiative)  stars,  the  slope  of  this  line  is  about 
0.25,  and  for  intermediate-mass  stars  (~1  Me,  mostly  convec¬ 
tive)  the  slope  is  about  0.06,  the  two  values  being  remarkably 
close  to  the  radiative  and  convective  exponents  cited  above. 
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Quantity  Value  Error 

Rotating  Model  Parameters 

Fa,  (km  s'1) .  273  13 

Leq.fi3  (km  s  ‘) .  374  3 

0  (day-1) .  2.71  0.11 

fi/fday-1) .  3.01  0.06 

Tp  (K) .  8740  140 

/,...  (K) .  6890  60 

Rp  (fl3) .  1.636  0.022 

/?cq  («0) .  1.988  0.009 

&mm  (mas) .  3.056  0.047 

9msa  (mas) .  3.598  0.017 

log  L  (L0) .  1.027  0.011 

log  gp  (cgs) .  4.266  0.012 

loggeq  (cgs) .  3.851  0.035 

Nonrotating  Parameters15 

M  (M0) .  1.791  0.018 

R(R&) .  1.652  0.022 

log  g  (cgs) .  4.256  0.002 

reff  (K) .  8200  98 

Xc .  0.607  0.019 


Note. — Uncertainties  due  to  the  parallax  have  not  been  included 
in  the  errors. 

a  Rotating  at  breakup  but  with  the  same  mass  and  polar  radius. 

15  The  parameters  of  a  nonrotating  star  from  the  Geneva  grid 
(Schaller  et  al.  1992)  which  would  reproduce  the  (corrected)  lumi¬ 
nosity  and  polar  radius;  see  text. 

Working  with  the  interiors  codes.  Claret  (2004)  is  able  to  eval¬ 
uate  this  exponent  at  each  point  in  the  evolutionary  paths  of  mod¬ 
els  covering  40  >  M@  >  0.08,  offering  the  results  as  appropriate 
exponents  to  use  in  rotating  stars  and  stars  in  close  binary  sys¬ 
tems  over  their  entire  evolutionary  lifetimes. 

This  is  a  constructive  suggestion  for  the  thorny  problem  of 
choosing  an  appropriate  gravity-darkening  law.  However,  we 
are  not  fully  convinced  of  the  leap  of  going  from  deriving  a  quan¬ 
tity  based  on  evolutionary  changes  to  using  it  to  describe  the  ef¬ 
fects  of  rotational  distortion.  One  might  make  the  case  for  small 
amounts  of  rotation  or  if  it  could  be  shown  that  rotational  distor¬ 
tions  and  evolutionary  effects  were  close  to  being  homologous 
transformations  from  one  to  the  other.  But  rotational  distortions 
are  not  homologous  to  evolutionary  changes,  and  it  is  not  at  all 
clear  how  well  these  “interpolations”  work. 

The  observations  reported  here  bear  somewhat  on  this  prob¬ 
lem.  We  have  tried  converging  our  Roche  models  using  the  value 
/ 3  =  0.09  adopted  by  Reiners  &  Royer  (2004).  The  results  are 
shown  in  column  (5)  of  Table  1.  The  V2  data  set  was  used  so 
the  comparison  is  with  column  (4).  To  achieve  the  degree  of 
asymmetry  found  in  the  triple  phase  data  with  this  low  exponent 
value,  the  rotation  parameter  is  forced  to  near  critical  rotation, 
to  ~  0.978.  In  turn,  the  predicted  projected  rotational  velocity 
v  sin  i  ~  295  km  s_1  conflicts  with  the  observed  value,  and  the 
predicted  color  is  significantly  redder  than  observed.  Further,  the 
reduced  y2  is  significantly  worse  for  this  fit. 

We  feel  it  is  premature  to  use  these  observations  to  derive  a 
“best”  value  of  the  gravity-darkening  parameter  until  the  re¬ 
maining  visibility  residuals  are  better  understood.  Since  it  is  the 
phase  measurements  that  we  trust,  which  are  sensitive  to  asym¬ 
metries  in  brightness  across  the  disk,  and  not  the  visibility  am¬ 
plitudes,  it  does  appear  that  the  von  Zeipel  (1924)  value  for  that 
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Fig.  4. — One  of  the  observed  V2's  for  the  AE-AW  baseline.  The  effect  of  the 
strong  H/3  feature  in  the  a  468.3  channel  is  clearly  evident  and  well  matched  by 
the  model. 


parameter  is  superior  to  the  value  adopted  by  Reiners  &  Royer 
(2004). 

Moreover,  as  we  describe  in  §  6.2,  our  results  indicate  that 
Altair  has  hardly  evolved  from  the  zero-age  main  sequence.  The 
Claret  (1998,  2004)  tables  give  values  for  f3  in  essential  agree¬ 
ment  with  the  von  Zeipel  (1924)  result  during  this  phase  of  evo¬ 
lution.  Thus,  Altair  does  not  provide  a  test  of  those  tables  but 
does  highlight  a  difficulty  in  applying  the  technique  proposed  by 
Reiners  &  Royer  (2004)  for  finding  total  velocities,  namely,  the 
difficulty  of  obtaining  a  priori  reliable  estimates  for  (3. 

5.4.  H/3 

In  Figure  4  we  show  the  blue  squared  visibilities  of  the 
AE-AW  baseline  for  the  H.A.  =  0b83  observation.  The  notable 
feature  at  486  nm  is  H (3;  the  agreement  with  the  calculations 
shown  in  this  scan,  and  the  others  not  shown,  is  striking.  This  fea¬ 
ture  is  nearly  centered  in  the  486.3  nm  channel.  In  contrast,  Ha, 
which  has  a  smaller  equivalent  width  and  is  split  between  chan¬ 
nels  at  2665.4  and  2648.7,  is  much  less  noticeable.  The  main 
reason  for  the  feature  being  reflected  in  the  amplitudes  is  the  re¬ 
duction  in  the  limb-darkening  coefficient,  the  star  appears  to  be 
more  like  a  uniform  disk  at  this  wavelength,  and  thus  the  visibility 
is  reduced.  This  close  agreement  is  a  nice  confirmation  of  the  de¬ 
tails  of  the  model  fits. 

6.  CONCLUSIONS 
6.1.  Imaging  Altair 

The  primary  result  of  these  observations  is,  we  believe,  the 
first  detection  of  asymmetric  surface  intensities  on  the  surface  of 
a  star  induced  by  rotation.  Although  we  have  imposed  a  model 
on  the  data  and  fit  the  model  parameters,  the  simple  conclusion, 
first  reported  by  Ohishi  et  al.  (2004),  is  that  the  surface  of  Altair 
displays  an  extremely  asymmetric  intensity  distribution  and  that 
the  asymmetry  is  consistent  with  that  expected  from  the  known 
high  rotation  and  with  the  previously  reported  oblateness  (van 
Belle  et  al.  2001).  We  have,  in  effect,  imaged  the  surface  of  an 
A  star. 

In  Table  2  we  summarize  various  physical  parameters  for  the 
adopted  model  (Table  1,  col.  [2]).  Most  quantities  should  be  self- 
explanatory.  Subscript  B  refers  to  the  model  if  it  were  rotating  at 
breakup.  The  angular  diameter  6m ;n  is  for  the  projected  minor 
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axis,  while  (?max  is  for  the  projected  major  axis,  i.e.,  the  angular 
diameter  of  the  equator. 

6.2.  The  Status  of  Altair 

Knowing  now  the  rotational  state  of  Altair,  we  can  better 
answer  fundamental  questions  such  as  its  evolutionary  status.  Over 
the  years,  Altair  has  been  classified  as  A7 1V-V  (see,  e.g.,  Johnson 
&  Morgan  1953),  the  luminosity  class  usually  indicating  an  object 
slightly  past  the  end  of  its  main-sequence  evolution,  while  the 
spectral  type  is  that  of  a  star  having  an  effective  temperature  in  the 
vicinity  of  7800  K  (see,  e.g.,  Erspamer  &  North  2003;  van  Belle 
et  al.  2001).  In  the  context  of  analyzing  Altair’s  pulsations,  it  is 
important  to  know  in  detail  its  evolutionary  state,  mainly  the  extent 
of  its  core,  and  to  remove  any  biases  that  might  be  introduced  by 
rotation. 

Fortunately,  in  the  context  of  rigid  rotation,  this  is  not  so 
difficult.  Early  results  (e.g.,  Sackmann  1970)  showed  that  two 
quantities  were  relatively  insensitive  to  the  effects  of  rotation: 
polar  radius  and  total  luminosity.  The  Roche  model  fits  give  po¬ 
lar  radius  directly,  while  it  is  a  straightforward  matter  to  calculate 
the  total  luminosity  (e.g.,  Collins  1963). 

These  two  quantities  are  not  perfectly  conserved.  In  the  range 
of  interest,  3  M0  >  M  >  1.4  Me,  both  quantities  decrease  with 
increasing  rotation,  approximately  in  proportion  to  u>2,  reaching 
a  maximum  correction  of  about  6%  in  luminosity  and  1.5% 
in  radius  (Sackmann  1970).  For  stars  in  the  neighborhood  of 
1.8  M0,  rotating  with  w2  ~  0.8  we  find  from  the  Sackmann 
(1970)  calculations  that  the  nonrotating  star  would  be  4%  more 
luminous  and  1%  larger  than  our  deduced  polar  radius. 

To  estimate  the  parameters  of  the  appropriate  nonrotating 
star  we  have  used  the  evolutionary  tables  by  the  Geneva  group 
(Schaller  et  al.  1 992).  These  models  were  calculated  with  modest 
convective  overshoot  (0.2 Hp).  We  have  used  the  grid  for  a  com¬ 
position  ofX  =  0.68,  Y  =  0.30,  Z  =  0.02. 

The  quantities  given  in  Table  2  as  “nonrotating”  are  those 
estimated  from  the  Schaller  et  al.  (1992)  models  and  are  quite 
striking.  The  last  entry  in  Table  2  is  the  mass  fraction  of  hydro¬ 
gen  remaining  in  the  core.  This  is  to  be  compared  to  a  starting 
value  ofXc  =  0.68.  Altair  is  almost  on  the  ZAMS. 

6.3.  Chromospheric  Indicators 

As  mentioned  in  §  1,  Altair  is  one  of  two  A7  objects  (the  other 
being  a  Cep)  in  which  certain  ultraviolet  emission  lines,  taken 
as  indicators  of  chromospheric  temperature  inversions,  are  seen. 
No  objects  of  earlier  spectral  type  are  known  to  show  these 
features,  and  it  is  usually  argued  that  these  therefore  represent 
the  hottest  photospheres  for  which  convection  is  still  capable  of 
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creating  such  temperature  profiles.  The  model  for  Altair  adopted 
here  calls  that  conclusion  into  question.  As  shown  in  Table  2, 
Altair  has  a  broad  swath  of 6900  K  gas  at  its  equator,  which  is  the 
likely  source  of  the  strong  convection.  We  also  note  the  recent 
announcement  (van  Belle  et  al.  2006)  that  substantial  oblateness 
and  gravity  darkening  have  been  found  in  a  Cep,  suggesting  sig¬ 
nificant  amounts  of  cool,  convective  gas  in  that  object  as  well. 

6.4.  A  6  Scuti  Star 

As  we  also  indicated  in  §  1 ,  rather  than  just  being  another 
rotating  A  star,  Altair  may  prove  to  be  a  valuable  laboratory  for 
examining  the  internal  rotation  state  of  a  star  with  a  predomi¬ 
nantly  radiative  envelope.  Buzasi  et  al.  (2005)  have  announced 
the  discovery  of  6  Scuti  pulsations  in  Altair  and  have  identified 
several  of  the  periods  with  frequencies  mostly  in  the  range  15— 
18  day-1.  Interestingly,  two  of  the  frequencies  reported  by  Buzasi 
et  al.  (2005),  3.526  and  2.57  day-1,  were  significantly  lower,  the 
latter  being  quite  close  to  the  rotational  frequency  we  have  derived 
(e.g.,  Table  2). 

This  was  immediately  recognized  as  providing  a  potential 
probe  of  the  interior  structure  and  particularly  the  rotation  law, 
and  an  attempt  has  been  made  to  identify  and  model  the  modes 
(Suarez  et  al.  2005).  Unfortunately,  not  knowing  the  rotational 
state  of  the  star  and  making  the  assumption  that  equator-on  was 
the  most  likely  orientation,  Suarez  et  al.  (2005)  adopted  a  total 
rotation  significantly  lower  than  what  now  seems  likely.  Other 
effects  of  this  choice  included  identifying  the  evolutionary  state 
as  being  substantially  more  advanced  than  we  believe  is  the  case. 
As  is  clear  from  their  results,  rotational  velocities  above  the  1 80— 
240  km  s-1  range  they  investigated  lead  to  large  changes  in  the 
oscillation  modes,  making  mode  identification  difficult.  It  will 
not  be  an  easy  task  to  tap  the  information  being  provided  by  Altair. 


The  referee,  Gerard  van  Belle,  asked  a  number  of  penetrat¬ 
ing  questions  that  significantly  improved  this  work.  Notably,  he 
wondered  about  the  effects  of  rotation  on  our  calibrators,  which 
led  to  the  discovery  that  we  could  measure  those  effects  in  (j  Aql. 
The  NPOI  facility  is  a  collaboration  between  the  Naval  Research 
Laboratory  and  the  US  Naval  Observatory  in  association  with 
Lowell  Observatory,  and  was  funded  by  the  Office  of  Naval 
Research  and  the  Oceanographer  of  the  Navy.  This  research  has 
made  use  of  the  SIMBAD  literature  database,  operated  at  CDS, 
Strasbourg,  France,  and  of  NASA’s  Astrophysics  Data  System. 

Facility:  NPOI 
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